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A NONLINEAR TIME COMPACTNESS RESULT 
AND APPLICATIONS TO DISCRETIZATION 
OF DEGENERATE PARABOLIC-ELLIPTIC PDFS 

B. ANDREIANOV, C. GANCES, AND A. MOUSSA 


Abstract. We propose a discrete functional analysis result suitable for prov¬ 
ing compactness in the framework of fully discrete approximations of strongly 
degenerate parabolic problems. It is based on the original exploitation of a 
result related to compensated compactness rather than on a classical estimate 
on the space and time translates in the spirit of Simon (Ann. Mat. Pura Appl. 
1987). Our approach allows to handle various numerical discretizations both 
in the space variables and in the time variable. In particular, we can cope 
quite easily with variable time steps and with multistep time differentiation 
methods like, e.g., the backward differentiation formula of order 2 (BDF2) 
scheme. We illustrate our approach by proving the convergence of a two-point 
flux Finite Volume in space and BDF2 in time approximation of the porous 
medium equation. 


1. Introduction 

There exists a large variety of numerical strategies for discretization of evolution 
PDEs. Proofs of convergence of many different numerical schemes often take the 
following standard itinerary (see, e.g., [27]). Given a PDE, discrete equations of 
the scheme are rewritten under a form reminiscent of the weak formulation of 
the continuous problem; stability estimates are obtained, which ensure bounds in 
appropriate (possibly discretized) functional spaces; eventually, sufficiently strong 
compactness arguments permit to pass to the limit in the discrete weak formulation. 

In many applications, including degenerate parabolic equations of various kinds, 
the question of strong compactness in spaces of sequences of approximate solu¬ 
tions is a cornerstone of such proofs. While “space compactness” is usually obtained 
by suitable a priori estimates of the discrete gradients involved in the equation, 
“time compactness” is often obtained by explicitly estimating or time trans¬ 
lates in the spirit of Alt and Luckhaus [Ij. This step is equation-dependent, and 
has to be reproduced for each problem. We refer to [27, 29] for the main ingredients 
of this already classical argument used in a number of subsequent works, to [5, 3] 
for some refinements, and to [15, §3] for a shortened version of the argument. 

The same question of time compactness often arises in existence analysis for 
PDEs, in the continuous framework. Along with the technique of [1], there exist 
several ready-to-use results that yield space-time precompactness of a sequence of 
(approximate) solutions {un)n- They are based on the two following ingredients: 

(A) estimates in sufficiently narrow Bochner spaces 

ensuring uniform in n bounds on space translates of the sequence (!/.„)„; 

{B) some very weak (in the space variable) estimates on the sequence {dtUn)n- 
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Then, different arguments permit to derive from [A) and {B) uniform in n estimates 
of time translates of (un)n and conclude that {un)n is compact in the appropriate 
space (e.g., as a consequence of the Frechet-Kolmogorov compactness criterion for 

spaces). This kind of result, in the abstract linear setting, is often called Aubin- 
Lions-Simon lemma [10, 37, 46], but the version we are interested in is also related 
to the early nonlinear version of the argument due to Dubinskii [24] (see also recent 
references [12, 20, 19]) and to the more recent formulation of Maitre [38]. Another 
related argument of nonlinear kind is due to Kruzhkov [36] . Further improvements 
were obtained by Amann in [2] for a refined scale of spaces (including Besov spaces 
for instance), and broached by Roubfcek in a rather general setting, see [44]. One 
observes that several closely related results co-exist, but the precise assumptions 
and conclusions of there results differ. Therefore, one can see the combination of 
properties {A)&l{B) as a “time compactness principle”, which can be made precise 
upon choosing a suitable functional framework or a suitable form of the estimates 
(A) and {B) (as one illustration, we refer to Emmrich and Thalhammer [25] where 
{B) is formulated as a fractional time derivative estimate). 

Further, discrete variants of different time compactness results have been already 
proved both for many concrete applications (mainly in the context of finite element 
or finite volume methods) and in abstract form: we refer in particular to Eymard et 
al. [27, 29] for Alt-Luckhaus kind technique for concrete applications, to Gallouet 
and Latche [31] for a discrete Aubin-Lions-Simon lemma, to Andreianov et al. [7, 4] 
for a discrete Kruzhkov lemma, and to Dreher and Jiingel [21] (see also [19]) for a 
discrete Dubinskii argument with uniform time stepping. 

The new result we are intended to present is based upon the technique of [39] 
which carries on rather easily to the discrete case. Our result is formulated as 
the ready-to-use Theorem 3.9 applicable to a large variety of numerical schemes 
(including variable time step and multi-step methods in time). It can be applied 
to a wide variety of strongly degenerate nonlinear parabolic equations. 

1.1. Compactness arguments for degenerate parabolic PDEs. In what fol¬ 
lows, fl is a bounded open subset of while T > 0 is an arbitrary finite time 
horizon. The cylinder D x (0,T) is denoted by Qt- 

Depending on the type of degeneracy of the underlying PDE, some of the afore¬ 
mentioned time compactness lemmas or techniques can be applied and the others 
fail to fit the structure of nonlinearities involved in the equation. To be specific, 
consider as the fundamental example the family of degenerate parabolic equations 

(1) = Az;-I-LOT, u G f3{v) 

where /3 C is a maximal monotone graph and LOT represent some lower-order 
terms, e.g., of convection and reaction kind. Definition and basic properties of 
monotone graphs are recalled in §1.2 for readers who are not familiar with this 
notion. The graph /3 can contain vertical and horizontal segments, which leads to 
elliptic-parabolic and (in presence of first-order convection terms LOT) to parabolic- 
hyperbolic strong degeneracy, respectively. In fact, most of the different time com¬ 
pactness arguments were developed for applications of the kind (1), with possible 
coupling to other equations. 

Let us assume for simplicity that we have a sequence of solutions to (1) with uni¬ 
form in n control of {vn}n in L^(0,T; iJQ(D)) and of {dtUn}n in L^(0, T; iJ“^(D)). 
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• Firstly, the Aubin-Lions lemma [10, 37] can be applied in this situation, pro¬ 
vided the graph /3 is bi-Lipschitz, i.e., if it is Lipschitz continuous with Lipschitz 
continuous inverse. This is the non-degenerate, uniformly parabolic case. Even in 
the power case j3{v) = sign(z;)|i;|“, the degeneracy at zero {a G (0,1), the porous 
medium equation) or the singularity at zero (a > 1, the fast diffusion equation) do 
not permit to apply the Aubin-Simon lemma. The classical generalization by Simon 
[46] of the Aubin-Lions lemma gives more precise compactness information under 
less restrictive estimates of {vn}n and of {dtUn}n, but it does not help to overcome 
the possible degeneracy of p. Degenerate cases require arguments of nonlinear kind. 

• The elliptic-parabolic degenerate case (i.e., the case where (3 is a continuous 
map) can be handled using Maitre’s lemma [38] or the Kruzhkov’s one [36, 9, 
3]. Application of each of these tools requires some additional assumptions such 
as the uniform continuity of (3 or the boundedness of ||u„||oo- The difficulty in 
application of Maitre’s lemma [38] consists in the choice of appropriate functional 
setting according to the behavior of the nonlinearity (3. It may require restrictions 
on the behavior of (3 at infinity and introduction of special functional spaces, e.g. 
of the Orlicz kind. 

• In what concerns discrete versions of the above general lemmas, the proof pre¬ 
sented in [36] is particularly simple to adapt to the discretized setting, indeed, it 
uses only the most natural norm for discrete solutions. Adaptation to the dis¬ 
crete setting of the linear compactness lemmas of [10, 37, 46] is presented in [31]; 
it requires the construction of the adequate discrete spaces and an ingenious re¬ 
formulation of the assumptions. It is feasible that also the nonlinear compactness 
lemma of [38] can be adapted to discretized setting using the approach of [31], but 
this question is beyond our scope. 

• Further, the elliptic-parabolic case (i.e., the case where /3 is a continuous map) and 
also the parabolic-hyperbolic case (i.e., the case where /3“^ is a continuous map) 
can be dealt with using the Alt-Luckhaus approach. It can be roughly described 
as estimating the integral over Qt of the product 

(2) (unix,t + t) - ^u„(x,t -I- r) - u„(x,t)^ 

by a uniformly vanishing, as r —>■ 0, modulus of continuity. Calculations leading to 
such estimate of (2) use the variational structure of the equations and the Fubini 
theorem; although they are standard, they appear to be equation-dependent and 
(in the discretized setting) scheme-dependent. 

Starting from (2) and estimates of Vu„, conclusions about relative compactness 
in L^{Qt) of {un}n (in the elliptic-parabolic case), respectively of {u„}„ (in the 
parabolic-hyperbolic case) can be derived, provided j3 (respectively, I3~^) is assumed 
to be Lipschitz continuous. Mere uniform continuity of j3 (respectively, j3~^) is 
enough for compactness, see [5]. 

Adaptation of estimates of the kind (2) to the discrete setting became a classical 
trend, starting from [27]. Yet the use of Fubini argument in the time-discretized 
setting brings lengthy technicalities, that are reproduced in an impressive num¬ 
ber of papers dealing with convergence of finite volume approximations to various 
parabolic problems. Therefore, our goal is to provide a black-box avoiding these 
computations. 
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1.2. Description of our approach in the continuous setting. To give an 
idea, in this paragraph we argue on equations (1). Let be (approximate) 

solutions of the problem. Given the structure of the equation, one can require that 
{vn)n is controlled in L^( 0 ,T;and {dtUn)n is controlled in the dual space. 
In the case /3 = Id, i.e. = Vn (with an immediate extension to bi-Lipschitz 
/3), this kind of assumptions is the basis of the Aubin-Lions-Simon kind lemmas; 
they are exploited for estimating space and time translates of the solutions, in 
order to apply the Frechet-Kolmogorov compactness argument. The idea of the 
compactness lemma we prove in this paper consists in justifying, under the same 
kind of assumptions, the relation 

(3) UnVn uv iu V'{Qt) 

(up to extraction of a subsequence), where u,v are the respective weak limits of 
Un,Vn (say, in L^{Qt))- Then we exploit this weak convergence property thanks to 
the monotonicity of (3. 

Motivated by application in porous media flows (see e.g. §1.3), we introduce a 
sequence (uJn)„ of L°° weights with L°° inverse converging almost everywhere to 
some limit uj. More precisely, we assume that there exists w, ZU > 0 such that 

(4) w < a;„(x) < uJ and ojn(x) —> oj(x) for a.e. x G fl. 

n—>00 

Let us first state a slightly modified version of the key technical tool [39, Prop. 1]. 
This result is reminiscent of the celebrated framework of Tartar-Murat (see [48, 40]), 
even though to our knowledge, there is no direct relation between the general theory 
of compensated compactness and the one of [39]. To simplify the statement, we 
exclude the cases g = 1 and <7 = 00 , where weak convergence should be replaced by 
the weak-* convergence (for {un)n, if (7 = 00 ; for {vn)n, if 9 = 1 )- 


Proposition 1.1. Let {uJn)n>i C L°°{n) be such that (4) holds. Let q € (l,oo) 
and p € [ 1 , 00 ), and let a > Assume that {un)n, {vn)n are two sequences of 

measurable functions on Qt such that (u„)„ is bounded in L‘‘{{0,T)]W^'°‘{fl)) and 
ivn)n is bounded in ((0,T);LP (f^))- Up to extraction of a subsequence, we can 
assume that and Vn weakly converge to u, v respectively in L'^((0,T); LP{fl)) and 
((0,T);LP (f^))- Assume that, in addition, (vn)n verifies the following uniform 
“weighted” weak estimate: 

(5) [ [ UJnVndtP < C'||Va;(p||cxj, Vtp G 'D{Qt). 

Jo Jfl 

Then 


( 6 ) 




0 JQ 


UJUV^p^ 


0 JQ 


ypGViQr). 


Observe the two main differences with [39, Prop. 1]. Firstly, weights ujn are intro¬ 
duced, which may be useful in PDFs modeling flows in heterogeneous porous media. 
For Proposition 1.1, the case of general a;„ follows from the particular one a;„ = 1 
(one can replace Vn by uJnVn), which is the one considered in [39]. Nevertheless the 
presence of general a;„ 7 ^ 1 will be more intricate to handle in the proof of Proposi¬ 
tion 1.4 (see below), this is why we keep track of uJn here. The second difference is 
the bound (5) which appears stronger than the one assumed in [39], in which (vn)n 
is only required to be bounded in some BV((0,T); space, where m can 
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be as large as needed. However, in our case we want to limit our considerations to 
m = 1 , because we will focus on numerical approximations and the information on 
higher-order in space discrete derivatives can be delicate to obtain. In this respect, 
the assumption (5) is the weakest assumption corresponding to m = 1. 

Before going further, let us recall the definition and a few basic properties of 
maximal monotone graphs. For more details, see for instance [16]. 

Definition 1.2. A monotone graph on R is a map P from R to the set 7^(]R) of 
all subsets of R such that (j/i — y2){xi — X2) > 0 for any (xi,X2) € and any 
( 2 / 1 , J/ 2 ) G P{xi) X P{x 2 )- It is said to be maximal monotone if it admits no non¬ 
trivial monotone extension. If P is a monotone graph, we denote by P~^ : R —>■ 7^(R) 
the map defined by y G P~^{x) x G P{y). 

Proposition 1.3. Given a monotone graph P, the following statements are equiv¬ 
alent: 

(i) P is a maximal monotone graph. 

(ii) P~^ is a maximal monotone graph. 

(Hi) for all A > 0, (Id-|-A/3)“^ : R —>■ R is a (single valued) non-expansive mapping. 

The “compensated compactness” feature ( 6 ) can be exploited in particular in 
the following context. 

Proposition 1.4. Let {un)n, (xn)n be two sequences of real valued functions on 
Qt weakly converging to limits u, v in L^{Qt). We assume that the limits u and 
V satisfy uv G L^(Qt). Assume that in addition, ( 6 ) holds. Let P be a maximal 
monotone graph with 0 G /3(0). If for all n, G P{un) a.e. in Qt, then v G P{u) 
a.e. on (0,T) x H. Moreover, up to the extraction of an unlabeled subsequence, for 
almost every (x, t) G Qt, either Un{:x-,t) —>■ u(x, t) orVn(x-,t) —>■ v(x.,t) and 

(i) Vn ^ V a.e. in Qt if P is single valued; 

(ii) Un ^ u a.e. in Qt if P~^ is single valued. 

For the general P such that neither /3, nor P~^ is single-valued one can describe 
the precise amount of strong convergence in terms of the support of the Young 
measures associated to the weakly convergent in L^{Qt) sequences {un)n, {vn)n- 
We defer to §2 the discussion of this general setting, the proof of Proposition 1.4 
and extensions of this result (see Remark 2.5). 

Remark 1.5. There are other ways to exploit the property (6). The way we propose 
in Proposition I .4 only brings L^ convergence, which is not optimal if additional 
properties of P are assumed (see, e.g., the exploitation proposed in [39, §4] or [22, 
§5] based on the convergence of some norm). However, the L^ convergence is the 
crucial fact, and it can be upgraded using equi-integrability bounds for stronger L^ 
norms. As a matter of fact, our result can be applied to a wide class of problems 
including strongly degenerate elliptic-parabolic and parabolic-hyperbolic problems. 

Proposition 1.4 is proved in §2, while §3 is devoted to the extension to the discrete 
setting of Proposition 1.1. A particular attention is paid in §3 to multistep discrete 
time differentiation like, e.g., Backward Differentiation Formulas (BDF). Finally, 
we apply our framework in §4 for proving the convergence of two-point flux in space 
and BDF2 in time Finite Volume approximation of the porous medium equation 

dtu — Au'^ = 0, q > 1. 
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But first, we give in §1.3 an example of how to use the combination of Proposi¬ 
tions 1.1 and 1.4 in the continuous setting in the case of the so-called Richards 
equation modeling the unsaturated flow of water within a porous medium. All the 
arguments we use can be transposed to the discrete setting, for example by using 
the numerical method proposed in [17]. 


1.3. A continuous example: Richards equation. We consider the Richards 
equation [43, 13] 

(7) uidts - div - pg)^ =0 in Qt, 

where the two unknowns, namely the saturation s G L°°{Qt', [0,1]) and the pressure 
head p, are linked by the capillary pressure relation s = S{p) where S' : R —>■ [0,1] 
is a nondecreasing function that satisfies S(p) = 0 if p < 0, S(p) > 0 if p > 0, and 
(1 — S) belongs to L^(R+). The porosity w G L°°{Qt) satisfies u<u}<uj a.e. in 
Qt- The intrinsic permeability field IK : fl —>■ A4£;(R) satisfies K(x) = K(x)^ and 

k\u\‘^ < K(x)m • u < Tcjuj^, Vm G R"^, for a.e. x G 


for some k, 7c > 0. The relative permeability k : [0,1] —7 [0,1] is increasing and 
satisfies fc(0) = 0, leading to a (weak) degeneracy at s = 0. The density p and 
the viscosity p are supposed to be constant, and g denotes the gravity vector. The 
equation (7) is complemented by the initial data 


(8) S|,^„=soGL-(r!;[0,l]), 
and the Dirichlet boundary condition 

(9) p|^,,, =pz5Gi7'(r!;R+). 

For discussions on more complex boundary conditions, see [45, 14]. 

In order to define properly the solution, we introduce the increasing Lipschitz 
continuous one-to-one mapping 


I R+ -7 R+ 

Ip HA /J’ yJk{S{a))<la, 


that is extended to the whole R as an odd function. 


Defiuitiou 1.6. A couple (s,p) is said to he a weak solution to (7)-(9) if ipip) ~ 
4 >{ pd ) G L^((0,r);i7Q(fl)), s = S{p) G (^([O, T]; Li(n)) with = sq, and 


( 10 ) 


ojdtS — div I ^ ^ ^ K ^V(^(p) — \/k{s)pg'j ) ~ ^ T^'{Qt)- 


We refer to [17] for a recent result of convergence of a carefully designed entropy- 
consistent nonlinear finite volume scheme for ( 10 ) in the discrete compactness 
framework developed in the sequel. Although our main interest is to describe a tool 
for numerical analysis of degenerate parabolic problems, here we limit ourselves to 
the continuous framework. Indeed, first, proofs of convergence of numerical schemes 
are very similar in their spirit to the proofs of stability of solutions with respect to 
perturbation of data, coefficients or non-linearities. Second, the results of Proposi¬ 
tions 1.1 and 1.4 are interesting already in the continuous framework. Therefore, 
here we focus on illustrating the structural stability feature (cf. e.g. [ 6 ]) of the 
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continuous problem (7)-(9). Let us show that a sequence (sn,Pn)n of solutions (in 
the sense of Definition 1.6) to equations 

(11) ojridts - div - pg)^ =0 in Qt, 

(playing the role of approximation of (7)) with the initial and boundary data (8)- 
(9) converges towards a solution of (7)-(9) in the sense of Definition 1.6. The only 
difference between (7) and (11) resides in introduction of approximate porosities 
(wra)ra>i that are supposed to satisfy properties (4). 

The derivation of the a priori bounds we propose here is formal. We refer to [18] 
for instance for a rigorous derivation on a closely related problem. 

Since s„ = S{pn) with S')®) = [0,1], we obtain directly that 0 < < 1 a.e. in 

Qt, ensuring the L^-weak convergence of an unlabeled subsequence (s„)„ towards 
some function s G L°°{Qt', [0,1]). 

Multiplying (formally) the equation (11) hy pn—po and integrating on Qt yields 
the estimate 



k{S{pn))\Vpn\'^dxdt 



\V(j){pn)\'^dxdt < C, 


Vn > 1, 


for some C independent on n. Thanks to Poincare’s inequality, one gets that 


\\4'{Pn)\\ L^{{0,T)-H^ (il)) — Vn > 1. 

In particular, there exists ^ G L^((0, T); with ^ — 4>{pd) G L^((0,T); TJg (fl)) 

such that (j){Pn) converges weakly in L^((0, T); i7^(D)) towards ^ as n —>■ oo. It 
follows from equation (10) that 

||wn9ts„||i2((o,T);ff-i(n)) ^ C, Vn > 1. 

Therefore, one can apply Proposition 1.1 to claim that 

(12) 0 JnSn4>iPn) in V'(Qt) as n —>■ oo. 

The sequences (srt)„ and (4>{Pn))n satisfy (('(pn) G l3(sn) for all n > 1, where 
P is the maximal monotone graph with single valued inverse = S o We 
can now apply Proposition 1.4. This first ensures that ^ G P(s), or equivalently, 
defining p : Qt —>■ K by p : that s = S(p). Second, the almost everywhere 

convergence of s„ towards s is ensured (up to a subsequence), so that one has 

Sn —> s and fc(s„) —> k(s) strongly in L’'((5t), Vr G [1,oo). 

n—^oo n—¥oo 

This is enough to pass to the limit in (11) and to claim that (s,p) is a solution 
to (7)-(9) in the sense of Definition 1.6. 


2. Proof of Proposition 1.4 

In order to avoid double integrals w.r.t. time and space, the functions we consider 
in this section are defined on an open O subset of . O will play the role that 
Qt played in the statement of Proposition 1.4. This allows in particular to write 
X G O rather that (x, t) G Qt- The weights oJn and oj are extended to the whole 
O, with 0 < w < Unix) < u and U!n{x) — >■ u{x) for almost all x G O. 
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Lemma 2.1. Let (un)^ and (vn)n sequences weakly converging in L^{0) 

towards u and v respectively Vn G (d{un) almost everywhere in O for all n > 1. We 
suppose that uv G L^(0), that UnVn € L^{0) for all n> 1, and that 


UJnUnVn^pdx 


ujuvipdx, 


to 


lo 


Then v S /3(m) almost everywhere in O. 


yip e v{Li). 


Proof. The proof relies on the fundamental property of maximal monotone graphs: 


(13) u S dom(/3) and v G /3(u) 

O {u — k){v — K)>0^ yk G dom(/3) and K G /3(fc). 

Let k G dom(/3) and K G /3(fc), then, for all p G T>{0) such that pPd, one has 


0 < 


Qr 


W„(Wn — k){Vn — itr)v?dx, Vn > 1. 


Due to the assumptions of the lemma, passing to the limit n —>■ oo one finds 


0 < 


Qr 


uj{u — k){v — K)pdTi., yp G VifD) wiih. p > Q. 


The above property is sufficient to claim that {u — k){v — K) > 0 a.e. in O. One 
concludes the proof of Lemma 2.1 thanks to (13). □ 


Keeping the notations of Lemma 2.1, we set Wn = Un + Vn- From Proposition 1.3 
we infer that A := (Id + and B (Id + /3)~^ are non-decreasing Lipschitz 

functions from K. to K satisfying furthermore A + B = Id, A(0) = 0 = i?(0) and 
one checks easily that Vn = A{wn) and = Biwn)- In the same way, we set 
w = V + u. From the assumptions of Lemma 2.1, we have Un ^ u and Vn ^ v 
weakly in L^{0), whence Wn ^ w for the same topology. Since = A(wn) and 
Un = B{wn), we may use the fundamental theorem on representation of weakly 
convergent in L^{0) sequences by Young measures (see [11, 33]), to deduce the 
existence of a Young measure {vx{-))xeO C Prob(R) (here Prob(IR) is the class of 
all probability measures on ffi.) such that for a.e. x G O there holds 

w{x) = / XdvxW, 

Jr 

and 

(14) v{x) = ( A{\)dvx{X), u{x)= f B{\)dvx{X). 

Jr Jr 

Lemma 2.2. Let (un)^ and (nrt)„ be two sequences as in Lemma 2.1. For almost 
all X G O, we have either A(A) = v(x) for v^— a.e. A G K or B{X) = u{x) for 
Vx— a.e. X G K. 

Proof. For G N, we denote by Ti : R+ —>■ K the truncation function defined by 
Te{r) = min{r, £}. Observe that Ti{A{wn)B{wn)) < UnVn due to the fact that 
AB > 0. The function w i—>■ Ti(A{w)B{w)) being continuous and bounded, we can 
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apply the fundamental theorem of [11] and claim that for all £ G N, one has 


/ (jjuvip da; = lim 
to ^ 


iOnUnVnV^X > llm 

, n—>oo 


U}nTe{UnVn)ipdx 


O 


= lim 

n—^oo 


u!nTi{A{wn)B{wn))(pdx = uj( Te{A{X)B{X))diyxW]‘pdx. 


lo 


lo 


Since this inequality holds for all £ G N, it also holds for the limit £—>■00 that we 
can identify thanks to the monotone convergence theorem. Bearing in mind the 
representation (14) of the functions u and v, this ensures that, for all Lp G ^{O) 
with </; > 0, there holds 

A(A)dPa;(A)'j ( f B(A)dz/a;(A)'j ipdx > f oj ( f yl(A)B(A)dp,r(A)'j tpdx, 

( / \t/]R / O / 

or equivalently (see, e.g., [34]) that 

ioifi J (^A{X) - A{fi)'^ (^B(X) - B{fi)^ dvx{X) dx < 0. 

The integrand of the above integral being nonnegative thanks to the monotonicity 
of A and B and the nonnegativity of Lp and w, it equals 0. Since Lp is arbitrary, we 
deduce that for a.e. x G O, one has 

- Mp)) {b{X) - B(^)) du^{X)dv^{p) = 0, 

whence (Al(A) — A{p,)){B{X) — B{fi)) = 0 for 0 i^x— a.e. (A, p.) G K^. Hence it is 
not difficult to conclude that either A is constant Vx— a-e., or B is is constant Vx— 
a.e., the value of the corresponding constant being fixed by (14). This concludes 
the proof of Lemma 2.2. □ 

Let us introduce the sets U and V defined by 

U ={x G O \ B = u{x) Vx — a.e.}, 

V ={x G O \ A = v{x) Vx — a.e.}. 

It follows from Lemma 2.2 that O \ (W U V) is negligible. 

Lemma 2.3. Under the assumptions of Lemma 2.1, one has 


u strongly in L^{U) and Vn 


V strongly in L^{V). 


Proof. We will prove that Un ^ u strongly in L^ifA), the proof of —>■ i; being 
similar. Since u G L^ifA), it can be approximated by simple functions: for all e > 0, 
there exist an integer and a simple function = X]i=i ^AEii where each Ki 
(1 < j < Jj) is a real value and each Ei is a measurable subset of U with |J- Ei = U, 
such that 

(15) I|m'- m||li(«) < e- 

Without loss of generality, we can assume that Ei is bounded if Ki ^ 0. For all 
i G {1,. ■. ,L}, the function w 1—>■ \B{w) — Ki\ is Lipschitz continuous, ensuring the 
uniform equi-integrability of the sequence {\B{wn) — Ri|)„, so that 


-u''|dx = y^ / \B{wn)-K^\dx —^ 

J^ n—^oo 


Ei 


|H(A) - Ki\dvx{X)dx. 
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Since Ei is a subset of U, we have 

/ \un — u'^\dx —^ / \u — u'^\dx, Ve > 0. 

Ju Ju 

It follows from the triangle inequality and from (15) that 

limsup / \un — u\dx < 2\\u'^ — u\\li(^h') < 2e, Ve > 0. 


n—¥c>o Ju 

Since e is arbitrary, we find that Un ^ u strongly in L^(W). □ 

The last lemma of this section focuses on the case where either /3 or is 
single-valued. 


Lemma 2.4. Under the assumptions of Lemma 2.1, and under the additional as¬ 
sumption that P is single-valued, Vn converges to v strongly in L^{0) as n ^ oo. 
Similarly, if is single valued, then Un converges to u strongly in L^{0) as 
n —>■ oo. 


Proof. Assume that /3 is single-valued, then the function B is (strictly) increasing, 
thus it is non-constant on any non-trivial interval. Therefore, Lemma 2.2 implies 
that, up to a negligible set, O = V, and one concludes by using Lemma 2.3. □ 

Remark 2.5. The assumption that maximal monotone graph /3 satisfy 0 G P(0) 
is easily dropped, indeed, it is enough to change in Proposition 1.4 the functions 
Un,Vn into Un — k,Vn — K respectively with {k,K) G p. It is also immediate to 
extend the result of Proposition 1.4 to a measurable in x family of maximal mono¬ 
tone graphs {P{x, ■))x^o (to fix the ideas, we can enforce measurability by requiring 
that the functions {x,z) >->■ {Id -\- P{x,-))~^{z), {x,z) ^ {IdP~^{x,-))~^{z) be 
Caratheodory), and 0 G P{x, 0) for a.e. x G O. At a price of some simple additional 
assumptions on {Pn)n, one can also consider the case of a sequence of convergent 
nonlinearities: Vn{-) G Pn{-,Un{-)). 

Finally, observe that the assumptions of weak convergence in O can be turned 
into LfJO) weak convergence assumptions; in this case, the conclusions in Lem¬ 
mas 2.3, 2.4 will turn into strong L]^^ convergences. 


3. Discrete “compensated compactness” result and its consequences. 

Our goal is now to derive a discrete counterpart of Proposition 1.1, namely 
Proposition 3.8. Combined with Proposition 1.4, it leads to Theorem 3.9 that can 
be used as a black-box. 

One needs to define discrete operators for defining a discrete function, its gradient 
and its time derivative. Rather than focusing on a particular numerical method, 
we concentrate on the fundamental properties a numerical method has to fulfill 
so that our result holds. This motivates the use of the so-called gradient scheme 
framework [23] for the spatial discretization in §3.1. Concerning the time dis¬ 
cretization, the approach discussed in §3.2 allows to consider either some one-step 
discretization methods, like for instance Euler and Runge-Kutta methods, or mul¬ 
tistep methods, like Backward Differentiation Formula (BDF). The main result of 
this section, namely Theorem 3.8, is stated and proved in §3.3. 

In what follows, LI is supposed to be a Lipschitz continuous open bounded sub¬ 
set of We restrict our attention to the case of cylindrical discretizations of 
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Qt := ri X (0,r), i.e., discretizations obtained thanks to a discretization of il and 
a discretization of (0,T). 

3.1. Spatial discretization. Concerning the space discretization, in order to be 
able to consider a wide range of possible numerical methods (including several 
Finite Elements with mass lumping and Finite Volume methods), we stick to the 
Gradient Schemes framework developed in [23]. 

Let m S N* be the number of degrees of freedom in space, we assume that there 
exist two linear operators 

TT^ : K™ ^ ^ 

such that the below assumptions (Axl)-(Ax3) hold. 


(Axl) For all p G [l,oo) and all m > 1, there exists a norm u i—>■ ||M||p,m on K™ 
such that the following assumption on the space translates holds: 


(16) 


lim sup sup -- = 0, 

ICI-S'O m>l 


Vp G [l,+oo), 


the function being extended by 0 outside of fl. 


In particular, a bounded w.r.t. the norm || • sequence 

a relatively compact sequence [TTmUm] in LP(Q). 

V / k>0 


c K™ yields 

V / k>o 


Remark 3.1. A classical choice for the norm || • ||p,m is 

||'^m|lm,p ■— ll'^m'^mllLP(Q) “t |1 ^m'^m |1 7 '^'^m G M , 

as suggested in [23]. Some spatial discretizations enjoy discrete Sobolev injections 
(see for instance the appendix of [28] for the case of Dirichlet boundary conditions 
and Appendix B of [7] for the Neumann ones), i.e., the property (16) holds if one 
sets 

||'^m|[m,p ■— ||Tm'^m||L5(0) “1” j] ^m'^m jj Lp (0)7^; VlZ-m G M 

for any q > Notice that q can be chosen strictly smaller than p in that case, 

and that on the contrary to the notation adopted in [23], the subscript p in jj • 
does not necessary refer to IF^’^(n), but it refers to the more general property of 
“uniform in m compactness” of the operator Tim ■ (M™, Jj • l|m,p) LP(fl). 


(Ax2) Let Um = (um,»)i<j<^ and Vm = be two vectors of R"* such 

that Vm G P{um), i.e., Vm,i G (i{um,i) for all i G {1,... ,m}, then Tr^rtm G 
)d{TrmVm) for a.e. x G fl. 

We should stress that this assumption is restrictive: typically, it is fulfilled for 
piecewise constant discrete solutions produced by reconstruction operator and 
therefore our analysis is suitable for finite volume methods and for the simplest 
finite element methods (in particular, methods with mass lumping). 

Before formulating the last assumption, we need more notation. In the sequel, 
{^m)m>i ^ L°°(r2) denotes a sequence for which there exists w, w > 0 such that 

(17) uj < ojm < w, a.e. in 11, Vm > 1. 
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We additionally assume that there exists w € such that 

(18) ujm —t w almost everywhere in fl. 

m—¥oo 

Obviously, one has u<uj<u} a.e. in O. 

Given ip G our last assumption is formulated in terms of the subset 

Vmi'-p) of defined by 


Vmi’f) 



^ m'^ m'^ r 


(p)dx = 0 , MUra 



(Ax3) There exists a linear operator Pm '■ C'^ (fl) —>• R"* such that Pm^ G 
and such that there exists C not depending on m such that 

\\VmPmp\\L^in) < Oilv<^|Uc=o(o), yp G c^m. 

The set Vm{p) is the preimage under tt^ of the L‘^^{n) projection of p on the 
range of It is not empty, and it reduces to a singleton in the particular case 
where Vm t \\T^m'Vm\\L'^{Q,) is a norm on R”^. In the latter case, the linearity of 
Pm is automatic since the range of tt^ is a finite-dimensional subspace of 
Since PmP G Tmip), there holds 

(19) [ UJm'^mUm{TrmP7 nP-p)dx = 0, G 0))° (fl), G R™. 

Jn 


Remark 3.2. Since we focus here only on compactness properties, we do not re¬ 
quire the reconstruction operators iTm and to fulfill the natural consistency and 
conformity relations that are required in [23] for proving the convergence of the 
methods based on gradient schemes. For instance, in problems that are posed in 
a Sobolev space W^'^, the consistency and conformity relations read, respectively: 
\/p e w'^'P{n), 

(20) Sm{p) ■■= min {W'lTm'Vm - p\\LP(n) + mVm - '^p\\Lp(Q.y) —S' 0; 

fmGR"* \ / m—)-oo 


and 'ip G C^(n; 


(21) Wm{p) ■■= max 


Vr, 


I m'^m ■ P + r„divV3)dx 
JQ 


0 . 


3.2. Time discretization. 


3.2.1. Extension of the spatial reconstruction operators. A time discretization of 
(0,T) consists in a subdivision 0 = to < ti <■■■< tn = T oi the interval [0,T]. 
For fc G {1,..., n}, we denote Atn,k = tk - 4-1 and A4 = maxi<fe<„ At^^k- 
This allows to extend the operators : R™ —>■ and : R"* —t 

(L°°(0))‘^ into 

; K-x("+i) ^ L°°{Qt) and V” : R™x(”+i) ^ {L°°{QT)f 

( \0<fe<n / \ 

m”’,) G R™x(n+1) ^ / n.M g 

R™ for all fc G {0,. .., n}. Then for fc G {1,..., n}, we set 
(22a) 7r>:;(.,t) and t) = 


if t G (4-1,4], 
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and (formally, since {0} is of zero measure in [0,T]) 

(22b) 7r>:;(.,0)=7rrf and V>”(.,0) = Vrf. 

Thanks to these definitions, we can define a semi-norm || - \\p^m,q,n on by 

/ n \ 1/9 

(23) IKII™.n= ■ 


Let us now extend the operator Pm to the case of time dependent functions. We 
introduce the linear operator ■ <^([0, T]; (^“(n)) —>■ defined by 


(24a) (P»" = P„v=(-,4-i), Vfc e n}, V(^eC([0,T];C“(f2)), 


(24b) 


(P^^f = Pm^{; 0), e C([0, T]; L^(L!)). 


It results from Assumption (Ax3) that there exists C > 0 depending neither on m 
nor on n such that 


(25) 




< CII Vy. G CT{Qt). 


3.2.2. The one-step discrete differentiation operator. It remains to define a recon¬ 
struction operator 5m in order to approximate the time derivative of the function 
. Since the case of the one-step time differentiation operator plays a funda¬ 
mental role in the analysis carried out in this paper, we first define 5m '■ —>• 

L°°{Qt) by 


(26) 5lu-m{;t) = 




,,A ;—1 


At 


n^k 


iftG (4-1,4], V< G 


The one-step time differentiation operator 6 m enjoys the following particular 
consistency property. 

Lemma 3.3. For all ip such that p{-,0) = (p{-,T) = 0, and 

for all Urn G one has 

[[ {‘^mdmKn'^hPliT + ^mT^mKadtp) dxdt = 0. 

J J Qt 

Proof : Thanks to the definition (22) of the reconstruction operator iTm, using the 
classical summation-by-parts procedure one has 




Qt 


'^u^dtpdyidt =Y^ 








n p 

= y' / UJmT^mU'^'" ipi-,tk) - pi-,tk-l))djC 

n p 

= -U^^)dK. 
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Thanks to (19), we obtain that 





fc=i 


Jn 


The result of Lemma 3.3 stems from the definition (24) of the operator P'^ and 
from the definition (26) of the operator 6 ^. □ 


3.2.3. From one-step to multi-step operators. The key idea of our reduction argu¬ 
ment proposed below is to represent multi-step differentiation operators as linear 
combinations of Euler backward differences with shifted time: for instance, 

- 2u'=-i -b = -(u'= - u'=-i) - - u'=-2). 

2 2 2 ^ ' 2 ' 

We introduce the appropriate matrix formalism for computations based on this 
idea; it will be exploited in Lemma 3.10 and in Proposition 4.5. Let us represent 
G ("■+!) by the matrix 


(27) 


/ 


= 


,n,0 


^m.l 


e M 


n+l,mV 


Then the operator 6 ^ corresponding to one-step discretization introduced in (26) 

where we have set 


can be rewritten as o ] 

/ 1 0 


(28) 


At„ 

0 


Atn 


Atn 


Atn 


0 - 


Atn 


0 \ 
0 


At. 


b-J 


G M 


n+H 


Remark 3.4. Since is only used on the negligible set {t = 0} of [0,T] in the 
reconstruction defined in (22), the choice of the first line of the matrix M„ 

is arbitrary. The particular choice we did in (28) leads to a matrix M„ that is lower 
triangular and invertible. 


Denoting by ID)„ and Tji the matrices of A4re+i(]R) defined by 


/I 0 . o\ 

-1 1 : 


: 1 0 
Vo ••• 0 -1 1/ 


/I 0 ••• 0 \ 

0 Ati ■ ■. : 

: 0 

yo • • • 0 At„/ 


we get that M„ = 















NONLINEAR TIME COMPACTNESS RESULT AND APPLICATIONS 


15 


In what follows, we restrict our study to the discrete time differentiation opera¬ 
tors (5” that can be defined by 


(29) 


Xn « n _ n n 


Vm” e R™x("+i) 


for some lower triangular invertible matrix M„ belonging to AIn+i(R)- Following 
the discussion of Remark 3.4, we can enforce the first line of M„ — denoted with 
subscript 0 in accordance with (27) — to be equal to the first line of M^, namely 


(30) 


(Mk) =1 , (Mn) = 0 if fc e {1 ,... ,n}. 
V /o,o V /o,fe 


Requiring that M„ is a lower triangular matrix means that for approximating dtu 


on (tfc_i, tfe], one can only use the vectors 


o<e<k 


C 


which is fairly natural. 


Requiring that M„ is invertible means that, knowing the initial value 0) = 

TTmiu^^) of TT^u^ and its approximate time derivative S'^u^, one can reconstruct 


We require a last very natural property on that is supposed to vanish 
on constant w.r.t. time vectors of discrete unknowns. More precisely, let = 




0 <k<n 


^mx(n+i) such that ^ for all k G {1 ,... ,n}, then 


S^u'^ = 0 a.e. in Qt- This amounts to assuming that 


(31) 


n 

E 

e=o 


(m„) 


= 0 , 


k,e 


Vfc G {1,..., n}. 


Lemma 3.5. Under the above assumptions, there exists a unique invertible lower 
triangular matrix A„ G AI„(M) such that 


(32) 


— T 

■- 71 


/I 0 • • • 0\ 
0 


V 0 

Proof. The fact that A„ is lower triangular and invertible follows directly from the 
fact that T„, M„ and are. The only thing to be checked is that the first column 
of An is equal to (1,0,..., 0)^. It is first easy to check that 

/I 0 ••• 0\ 


ft-i 


1 


Vi 


0 

1 1 / 


The property ^A„^ = 1 follows from the particular choice (30) of the first line 


of M„, while the property (A„ ) = 0 for fc > 1 follows from (31). 

V / k,0 


□ 


With the above formalism, we thus focus on discrete time-differentiation opera¬ 
tors of the form 


(33) 




Vm” G R™x(n+1)^ 
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with An of the form (32). Note that the one-step differentiation enters this frame¬ 
work. In this case, the matrix A„ reduces to the identity. 

Remark 3.6. In the particular case of a uniform time discretization, i.e., Atn^k = 
T/N = Atn for all k € {1,..., n}, it results from the particular structure of the 
matrix An (see Lemma 3.5) that An and T„ commute, so that M„ = A„M„. 


The assumption we make on the discrete time-differentiation operator is: 

(At) Let 6 (!n : —>■ be an operator of the form (33). where A„ S 

AlTt+i(M) and A„ G AIr!,(]R) are the matrix defined by (32). We assume 
that there exists C not depending on n such that ||A“^||^ < C, where || • ||i 
denotes the usual matrix 1-norm, i.e., ||]B||]^ = maxi<j<„ ELi for all 
1 G 7W„(R). 

Let us stress that this assumption is fulfilled by the one-step differentiation opera¬ 
tors with arbitrary time steps since A„ reduces to Identity. However, for multi-step 
methods like BDF2 described below, one may need to constraint the ratio between 
adjacent time intervals in order to guarantee that (Aj) holds. 


Remark 3.7. Let us illustrate (in the particular case of a uniform discretization) 
how to determine the matrix A„ prescribed by (32) for the so-called BDF2 scheme 
(see e.g. [47]^. The principle a such a method consists in an initialization with a 
one-step differentiation 




*m _ 

Atn 


*/tG(0,At„), Vm” G 


while is defined on the following time steps by: G Vfc G {2,..., n}, 

^^,n,k _ 9^,n,fc—1 i 2 \ 

2 ^ ' o ^m. \ 




2 rn 


Atn 


if t £ {{k — l)At„, kAtn). 


This leads to 

(x,t)l(o,At„)(t) 


and thus, in view of Remark 3.6, to 

/ 1 0 


+ (X,f- At„)^ l(At„,T)(i), 

• ••• 0 \ 


(34) 


1 3 

2 2 


It is easy to verify that 

ll(A„)-^| 


V 0 


= ; 1 - 


- 0 

0 -I V 

" 2 2 / 


< 


2 ’ 


Vn > 1, 


so that the BDF2 method with uniform stepping satisfies Assumption (At). 
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3.3. Discrete “compensated compactness” feature. We can now state the 
main result of this section, which is a discrete version of Proposition 1.1. 


Proposition 3.8. Let ^>1 (^m)m n>i discrete reeonstruction opera¬ 

tors as defined in ^S.2.1. We assume that Assumptions (Axl) and (Ax3) hold. Let 
[Slf] be a family of discrete time-differentiation operators of the form (33) 

satisfying Assumption {At). Let and (Um)m,rt>i be two families of vec¬ 
tors such that G ]^™x("-+i)_ {u!m)m>i C L°°{i}) be a sequence of func¬ 

tions such that (17) and (18) hold. We assume that there exist p and q in (l,c») 
and C depending neither on m nor on n sueh that, for all m,n> 1, 
a. < C, in particular n>i admits a weak limit (up to a 

subsequence) u in L‘‘{{0,T); LP{fl)) as m,n ^ oo; 
b- that )^ „>^ admits a weak limit (up to 

a subsequence) v in L‘‘ ((0, T); L^ (fl)) as m,n ^ oo; 
c. for all ipif, G ]R'"x(n+i)^ 

(35) ff 

J J Qt 

Then, up to a subsequence, one has 



V^dxdt 



LOUVip&X.dt, 


yip G V{Qt). 


Combining Proposition 3.8 with Proposition 1.4 we get the following result, 
which is the main result of our paper. 


Theorem 3.9. Keeping the assumptions of Proposition 3.8, we additionally sup¬ 
pose that Assumption (Ax2) holds, and that there exists a maximal monotone graph 
(3 such that G for all m,n>l. Then the weak limits u,v of 

and {T^mVfrfim n satisfy u G P{v) for a.e. (x,t) G Qt- Moreover, up to an unlabeled 
subsequence, 

(i) If (3 is single-valued, then, v a.e. in Qt as m,n ^ oo. 

(ii) If is single-valued, then —>■ u a.e. in Qt as m,n^ oo. 


With Propositions 3.8 and 1.4 at hand, the proof of Theorem 3.9 is straightfor¬ 
ward. Indeed, Assumption (Ax2) ensures that 


u- G /I(n(^) (i.e. G /3{u))^,}) 


tt” u” G ^(tt” m” ) a.e. in Qt, 


so that one can directly use Proposition 1.4. The remaining of this section will be 
devoted to the proof of Proposition 3.8 


3.3.1. Reduction of the problem to a semidiscrete situation. In order to prove Propo¬ 
sition 3.8, our strategy consists in reducing Estimate (35) into a semi-discrete esti¬ 
mate that will be easier to handle. This is the purpose of Lemma 3.10 stated and 
proved below. 

Lemma 3.10. Let be a time-differentiation operator as introduced in %3.2.3. As¬ 

sume that (A.x3) and (At) are fulfilled, and that (35) holds for all G 
then 

(36) [[ w^7r”M”dt(pdxdt < C'||V(p||ioo(Q,^), y^p €V{Qt). 

J J Qt 
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Proof : Let G K™x("-+i) be arbitrary, then define = (A„^) ip^, so that, 
for all fc G {1,..., n}, one has 


(37) = 


E (*»')« 




1=1 




L“(n) 


< 


m IIL°°(Q) — ll'^n WlW^ mPmW (Qt) ' 


The link (33) between 6 ^ and provides that 


ff ^)),dxdt = //' 7r”^”dxdt 

JJqt JJQt ^ ' 


'n 


k=l 


' Qt 


i At n 1 

-TT, 


At. 


At„ 


7rmv)))^dxdt 




^ I ^m'^m 

Q. 


ri,n,i _ «,n,^—1 
“'m 


At 


n,i 


J2 (^n)Gfc^"‘^m''dxdt. 




Since = A^^” , one obtains that 

m” dxdt = [[ w™(5” m” tt” dxdt, 

' ^ Qt 

which ensures together with (35), (37) and (At) that 

(38) ff a;„C«<^:;clxdt<C||V(),v^"||^,.(Q^). 

Qt 

Let p G 'D[Qt) be arbitrary, and set — —P^p so that, thanks to Assump¬ 
tion (Ax3), one has 

< C||V<p|Uo.(Q,), Vm,n > 1. 

Thus it follows from Lemma 3.3 and (38) that 

ff UJmTT^U^dtpdxdt = ff UJmSl^u'f^TT^p^dxdt <C\\Vp\\l'^(^q^), 

JJO'T JJot 


concluding the proof of Lemma 3.10. 


□ 


3.3.2. Proof of Proposition 3.8. The proof mimics the one of [39, Lemma 3.1]. Let 
p G C^{Qt)^ we denote a = dist(suppv?; 911). Let p G C')?°(R‘^) be such that 
p(—x) = p(x) > 0, for all x G K"^, such that suppp C ^^(0,1) and such that 
p(x)dx = 1. For £gN, £>l/a and x G K"^, one defines p^(x) = £'^p(£x), so 
that supppf C 5^(0, l/£), fgd pe{^) = 1 and p* pt G Cf°{QT), where * is the usual 
convolution w.r.t. the space variable x. 

One splits 

(39) f f {UJUV - UlmTT^ul^TT^vl^) pdxdt 
J J Qt 

= Ri{tj + R 2 {^, m, n) + i? 3 (^, m, n) + i? 4 (^, m, n), 
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where 


— [uujv — u{{ujv) * / 9 ^)](/?dxdi, 

J JQt 

R2i£, m,n)= [uHujv) * pi) - tt” m” i;” * p£)](pdxdt, 

J JQt 

Rsii, m,n)= ii” * pe) - m” tt” ii” ) * pe]pdxdt, 

J J Qt 

Ri{£, m,n)= [[ it” tt” t>” ) * pe - m” tt” t>” )](^dxdi. 

J JO^ 


^Qi 

Clearly, {u)v)^pz tends weakly to ujv in ((0, T); (11)) as i tends to +oo, leading 

to 


(40) 


lim \Ri{i)\ = 0. 

£—¥oo 


Since pi is an even function, using the Holder inequality and the bounds assumed 
in the statement of the proposition we find 


|i? 4 (£, TO, n)| = 


// * pi - p)dxdt 

J JQt 


(41) 


C, 


<c\\p *pi- </5||L~(Qr) < — II Vv?||i=o(Q^). 


In particular, i? 4 (f,TO,n) tends to 0 uniformly with respect to to and n as £ tends 
towards +oo. 

Concerning the term R 3 {£,m,n), one has for almost all (x, t) S Qt ■ 

Sm.nA^^t) ■= 7rmMm(x,t)[(w„7r”u” ) *p^](x,t) - [(w^TT^ m” Tt” t>” ) * Pf] (x, t) 


< 


t) - n” u" (x - y, t) 


Wm(x - y)n" n)),(x - y, t)pKy)dy- 


JBd,io,i/e) 

We deduce from Fubini’s theorem that 

|i?3(^,TO,n)| < ||<p||l“(Qj,) ||<S'm,rt,£||^l(Q^) 

< wIIpIIloo(q^) IItt” v" f \\TT^u^-TT^u^{--y,-)\\Li(LP)Piiy)dy. 

It results from Assumption (16) and from the assumption ||Mm||p,m,( 3 ,n < C that 

/ lkm«m (• -y,-)l|L«(Lr)Pf(y)dy ^ 0 

dB4(0,£) 

uniformly w.r.t to and n, so that 


(42) 


R 3 {£,m,n) —> 0 uniformly w.r.t to and n. 

^—>• + 00 


Let us now focus on controlling the term R 2 {£,m^n). Fix £>1, and define the 
family ^ C L^{Qt) by 


Z^m,n,£ ■— (.^'TTi'^rn'^rn) * 


Vto, n > 1. 
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The sequence (uJmTr^vl^)m,n is uniformly bounded in L^(Qt), then it follows from 
the regularity of pe that there exists Q > 0 (possibly depending on £, but neither 
on m nor on n) such that 

(43) < Q, Vm,n>l. 

On the other hand, given ip € C^{Qt), we have 


{dtZm,n,i,i’)v'iQT)MQT) = - // Zm,n,^l9t'0dxdf 

Iqt 


[[ * pt)dyidt. 

J JQt 


Therefore, it follows from Lemma 3.10 that 


(44) \{dtZm,nJ,'>P)v'{QT),V{QT)\ ^ C\\V {ip * pt)\\oa < Cl\\lp\\c([ 0 ,T]-L^{Q.)) 

for some Ci possibly depending on t but neither on m nor on n. As a consequence, 
for fixed £, the sequence {dtZm,n,i)jn n bounded in the space of finite Radon mea¬ 
sures on Qt- Along with (43) this ensures that the family {zm,n,i)^ „ is bounded 
in the space BV(Qt), thus yielding 

(45) izrn,n,e)m,n>i i® relatively compact in L^{Qt). 

Since Zm,n,e is piecewise constant in time and smooth in space, the map 


hm,n,£ ■ 




X l-> 


/(0,T) 


\dtZm,n,£{^i t)| 


is continuous on ft. Let Xm,n,e G be such that 

bm,n,£(^m,n,^) ^ 7> SUp ldn,n,^' 

^ xen 

Let ip G C^{{0,T)), and let fc G N with k > 1/d{^rn,n,i,d£V), then choosing ip : 
(x, t) !->• p{t)pkpx. — Xm,n,£) in (44) and letting k tend to -boo yields 

bm,n,^(Xm,n,£) ^ CV: VtU, 72 ^ 1, 

whence 

(46) \\Vm,nA\L°°in) < Ci, Vm, n>l, 

where Ci depends on £ but neither on m nor on n. In addition, the family 
{oJmT^^v^)^ ^ being bounded in L^{Qt) = (fl; L^((0, T))), the definition of 
Zm,n,i ensures that 

\\z^rn,n,l\\c(Q.\L'^((£),T))) — ^171,71 > 1. 

Therefore, for all x G fl, there exists t^rnni ^ (bj T) (possibly depending on x, to, 
n and £) such that 

(47) \zm,nA^^tAm,n,e)\ ^ 'im,n> 1. 

For all X € and almost all t € (0, T), one has 

t)\<\ A^^x,Tn,n,f.) \ 4 ” An,n,l{A — 

thanks to (46) and (47), whence 

(48) llZm^nAlL-^-ipT) < ^171,11 > 1 
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for some Ct, depending neither on m nor on n. Due to (48) and (45), we de¬ 
duce that the family is relatively compact in for all r € 

[l,-|-oo). In addition, (up to an unlabeled subsequence) Zm,n,e converges weakly in 
((0,T);LP (fl)), as m,n —>• oo, to the limit {ujv) * Choosing r = max{p',( 7 '} 
we see that, up to an unlabeled subsequence, Zm,n,i converges towards (luv) * pi 
strongly in ((0,T);LP (D)) as m and n tend towards -l-oo. Therefore, since 
tt” m” converges weakly in L'^((0,T); L^’(fl)) towards u, we can claim that 

(49) lim |i? 2 (^, 171,71)1 = 0, V£ > 1. 

m,n—^■oo 

Let e > 0, then using (40), (42) and (41), there exists > 1 such that 

|i?l(4)| + 1-^3(4,771,71)1 -f |i?4(4,771,7l)| < €, VtTI, 71 > 1. 

It follows from (39) and (49) that 


lim sup 

m,n—¥oo 



(, , „n ^.n ..n 


concluding the proof of Proposition 3.8. 


ojuv) (pdxdt 


< e, 


Ve > 0, 


□ 


4. Application to Porous Medium equation 

The goal of this section is to provide a new convergence result for a two-point 
flux in space (cf. [27]) and BDF2 in time (see Remark 3.7) approximation of the 
solution to the porous medium equation as an application of Theorem 3.9. More 
precisely, let D be a polygonal subset of with outward normal n, let T be a finite 
time horizon, let q > 1, then given uq € L^(D), we aim to approximate the solution 
u of 

{ 9(71 — A(|7t|^“^7t) = 0 in Qt = ^4 X (0,T), 

V(|7t|«-iM) • n = 0 on 9D X (0, T), 

u{-, 0) = uq on D. 

Such an equation has been widely studied in the last decades, see in particular [42, 
49]. In particular, it is well-known that the problem (50) admits a unique solution 
u G L°°((0, T); L^(D)) and such that g L^((0, T); iJ^(D)) (see e.g. [1, 41]). 

We propose a formally second-order accurate in both time and space Finite Volume 

scheme, and show the convergence of the corresponding family approximate solution 
towards the unique solution u as the discretization parameter tend to 0. 

4.1. Discretization of Qt- We require the spatial mesh to fulfill the so-called oth- 
ogonality condition, sticking to the definition of [27, Definition 3.1] for an admissible 
discretization {T,£, {:x.k)kgt) of D. More precisely, the domain ft is supposed to 
be split in a tessellation T of open polygonal convex subsets {K}k&t, such that 
[JkgT ^ ~ ^ iFnL = 0ifiF ^ L, where {K, L) G T^. Each control volume 
K gT is endowed of a so-called center xk G K^, and of interfaces (edges, if d = 2; 
faces, if d = 3) CT C dK contained in hyperplanes of the (d — l)-dimensional 

Lebesgue measure of an interface ct being supposed to be strictly positive. The 
intersection of the closure of two elements K,L gT is either an interface (denoted 


^This assumption is made in order to lighten the presentation; it is a classical issue to relax this 
assumption, requiring only that point from K to L (this is the case, e.g., under the Delaunay 

condition for simplicial meshes) and that the mesh size h be defined by hj- = maxxeT^i^^(^ ^ 
{xic}); we refer, e.g., to [8] for details. 
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by (Jkl), or a manifold of dimension less than d — 1, or it is empty; we define the 
set £ of the interfaces hy £ = {axL \ K,L G T}. For all K £ T, we denote by 
A/*: C T the set of the neighboring cells of K, defined by L G Afx iff ctrl G £. We 
assume that for all itf G T and L G Mr, the segment [xr,xl] crosses the interface 
aRL orthogonally. We denote by ukl = normal vector to aRL 

outward with respect to K and inward with respect to L. Note that the above 
definition of £ does not contain the interfaces lying on the boundary d£l of £1. We 
define by £ext the set of such boundary interfaces, and £R^ext the subset of £ext 
made of the interfaces of d£l. In what follows, we denote by ttir the d-dimensional 
Lebesgue measure of df G T, by uirl the (d — l)-dimensional Lebesgue measure 
of aRL G £, and by rrin the (d — 1) dimensional measure of ct G i^ext- We introduce 
the size hq- and the regularity pj- of the mesh T by setting 

/mKilxK-Xil diam(A:) \ 

51 fiT = maxdiam A , pt = max >-h -r 

ReT IXK-Xiiy 

Concerning the time discretization of (0,T), in order to simplify the presentation 

we restrict our study to the simple case of uniform time discretizations. Given 
n G N*, we denote At = T/n and tk = kAt for fc G {0,... ,n}. 

Given Uj- = {uk)r^t G we denote by tt-j-Uj- the piecewise constant 

function defined almost everywhere in £l by 

(52) TTq-Uq-{x.) = Ur if X G AT. 

For ulj- = (u|^)R:gr,fce{o,...,n} G ]R#'Gx("+i)^ and fc G {0,..., n}, we denote by 

= {u'k)ReT G 

We denote by the piecewise constant function defined almost everywhere in 

Qt by 

= Trq-u^’^{x) = u%; if X G df and t G {tk-i,tk]- 

In order to reconstruct a discrete gradient, we introduce the so-called diamond 
cells Dkl for url G £, which are open subsets of £l defined as the convex hull of 
xr,xl and aRL- Given vj- = {vr)r^j- G we define 

(53) Vrvr{x) = d^ —^n^L, Wx £ Dkl, ^aRL £ £■ 

\^R - XL I 

Using the geometrical identity 

^ /" j toxl|xr-xl| ^ 

(54) ineas{DKL) = dx =- - -, vsrl £ £, 

-Idkl “ 

we obtain that 

(55) l|Vr'urlli2(Q)d = d ^ TRL{vR-VLf, 

(TKLeS 

where trl = , ■ This leads to the following definition of the discrete norm 

II • ||x,r: Vur = {vr)r^j- g 

Il^rll2,r =IKr^^rlli2(Q) + || Vr'wrlli2(Q)d 

^'^{vR^mR + d ^ tkl{vr-vlY- 

ReT cTKLeS 
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We also define ^ by 


4.2. The Finite Volume scheme. The initial data uq is discretized into = 
^ where 

(56) f Mo(x)dx, yK G T. 

rriK Jk 

It follows from Jensen’s inequality that 

(57) IKrMrllL=(n) = < ||wo|li2(n). 


In the sequel, we denote 


iIj : 


KeT 


u I—>■ ImI'^ ^u. 


We use the implicit Euler scheme for determining = {u]^) i-C. 


(58) 


'^K 

At 


rriK + ^ tkl (^(mk) - iPiuD) =0, VX G T. 




As soon as 2 < fc < n, we use the so-called BDF2 scheme for determining = 
{^k)kgT^ that is required to fulfill: VAT G T, Vfc G {2,..., n}, 


(59) 


_ 27/^ ^ -I- —'ji^ ^ 

2 K K ^ (V’(mI-) - '0(wi)) = 0. 


At 


4.3. Main a priori estimates and existence of a discrete solution. Define 
the function 0 : R —)■ R by 


(60) 


pU _ 

(/)(u) = / ^J^'{a)da = 
Jo 


^ |u| ”*2 u, Vu G R, 


0 9 + 1 

then the Cauchy-Schwarz inequality yields 
(61) (a - b){ilj{a) - tp{b)) > {(j){a) - 4>{b)f , V(a, b) G 

Therefore, multiplying (58) by Atu]^ and summing over K € T, using 

(a-6)a>y-y, V(a,6)GR^ 

and the classical summation-by-parts, we find 


( 62 ) i {^kY'^k + At ^ tkl {(I){uk) - (t>{.ul)Y {urY^k- 

^ KeT CTKL ^ kgT 

In order to obtain an estimate on the following time steps, we use the inequality 

-26-1- a > ^ (a^ + (2a — b)^ — 6^ — (26 — c)^) , V(a, b, c) G M^. 


/ 1 
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Multiplying (59) by Atu’^ and summing over K G T and k G {2,... for some 
£ G {2,..n}, we get 

((^1-)^ + mx + J^AtJ^ tkl - (l^iupp 

KeT k=2 (Jkl 

* KeT 

Combining (62) and (63) and using (57) we find that for all £ G {1,..., n}, 

4 {'^kY'^k + {(I>{uk) - < 2||uo||i2(n), 

KeT k=l (Tkl 

leading to the following statement. 

Proposition 4.1. Let uVj- be a solution of the scheme (58)-(59), then there exists 
C depending only on uq and d (but not on the discretization) such that 

ll^r“rllL“((0,T);L2(n)) + \\'^T^i'^T)\\ 

In order to approximate the solution to the scheme (56), (58) and (59), one can 
use the iterative algorithm based on monotonicity proposed in [27, Remark 4.9], 
that converges towards a solution to the scheme. Therefore, there exists at least 
one solution to the scheme. The uniqueness of the discrete solution follows from 
a classical monotonicity property of two-point flux approximation (see e.g. [26]), 
leading to the following statement. 

Proposition 4.2. Let (T, £, (x-k)kgt) be an admissible discretization of LI, then 
there exists a unique solution to the scheme (56), (58) and (59). 

We need another a priori estimate before applying Theorem 3.9. This is the 
purpose of the following statement. 

Lemma 4.3. Let Uj- be the unique solution to the scheme (56), (58) and (59), 
then there exists C depending on LI, T, q, pj- and d such that 

lkrV'(“r)llLi(QT) + ll^rV'(“r)llLi((3T)‘^ — 

Proof. First of all, by adapting to the discrete framework the technical lemma [35, 
Lemma A.l], we can claim that there exists C depending only on uq, q, T and d 
such that 

(65) ll^r'/'(“r)llL2(Qi,) ^ 

Indeed, in the continuous case this property relies on the Poincare-Wirtinger in¬ 
equality, for which there exist discrete counterparts [32, 7]. It follows from the 
definition (60) of the function (j) that there exists C depending only on uq, q, T and 
d such that 




Li+HQt) 


< C. 


( 66 ) 
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The definition of the discrete gradient and the geometrical identity (54) 
ensure that 

n 

(67) |lV^V'(«r)llLHQT)=E^^ E ^al|^(«^)-V'(4)| 

k — 1 
n 

= E^^ E ^KLri\L\^{u)i) - 

k —1 <yKL^S 


where, for all fc G {1,..., N} and all cjkl G we have set 


if u\ ^ u'= 


L’ 


ifn)^ = <. 


Vkl = { 

>'( wa ) = VqIukI"^ 

In particular, the mean-value theorem yields: € £,yk > 1, 

0 < <V^{\u’k\^ + \ 

Using the Cauchy-Schwarz inequality in (67) we get 


||V;fl/-(M^)||i,(Q^) < |1V^<).«)|U.(Q,) ^ mKL{v'kL)"\^K-^L\ 

\fc=l (yKL^S 
n 

<cj2^t rriKL {\u'k\'^~^ + \ul\‘^-^) \XK - xl\ 

k — 1 CTKL^S 

< At ^ I E Wifilxif- xlI 

fe=i KeT XLeMK 

Using the regularity pj- of the mesh defined by (51), one obtains that there exists 
C depending only on uq, T, Q, q, d, and pj- such that 

Using (66), one obtains that 

\\^Ti^i^T)\\LHQ.Y<C. 

In order to conclude the proof of Lemma 4.3, it only remains to use again the 
discrete counterpart of [35, Lemma A.l] to obtain that ||'T^V'(“r)llLi(QT) — 

4.4. Compactness of the solution. This is the point where the main results of 
this paper are exploited. Let {%m£rm {^k)K eTm)m>i be a sequence of admissible 
discretizations of Qt such that there exists p* > 0 satisfying 

(68) lim hr^ = 0, sup pr™ < P*- 

m^oo 

For the ease of reading, we denote by G instead of the discrete 

solution to the scheme (56), (58) and (59) corresponding to the mesh TTn and the 
time step At„ = T/n. Similarly, we replace the notations 717 -^, Vt^) and 
by 7 r„, tt” , Vm and respectively. 

For all m > 1, the functions 

Um^hmUmWL^in) and u m ||'^m||2,m ■— II||+ II || 
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define norms on Let us check the assumptions of §3.1 on the space dis¬ 

cretization. It follows from the so-called space-translate estimate [27, Lemma 3.3] 
that 

+ 4) ~ '^TnVm\\L^{Q) 

lim sup sup --li-= (J, 

SO that (Axl) holds. Since tt^ consists in the piecewise constant reconstruction, 
Assumption (Ax2) also holds. Finally, observe that Um >■ \\T^mUm\\ 1 ^ 2 defines a 
euclidian norm on . It is easily checked that the linear operator Pm defined 

for (f e C^{^) by ^ with 

<PK = — / <p(x)dx, VA e Tm, 

rriK Jk 

is suitable in Assumption (Ax3). Indeed, the discrete gradient of Pm'P is given by 

T 7 r, ,.AA - ^L){y^K -y^L) w_ 

V mPm‘-P\p^) — a I , Vx € PkL^ VfT/CE € 

|xk - Xij^ 

Thanks to the continuity of ip, there exists {^K)KGTm such that S K and 
(px = V’(xr:) for all A G 7^. Therefore, it is easy to verify that 

\\VmPm^\\L<^(^Q) < d{l + 2p*)|| Vvj||ioo(o), Vm > 1, 

so that Assumption (Ax 3) holds true. 

Let p G C^{fl X [0, T)), then multiplying the scheme (58) (resp. (59)) by Atp'^^ 
(resp. l^tp'x^ for k >2) where 

= / </5(x,tfc_i)dx, VA G 7)n, Vn G {1,... ,n}, 

JK 

and summing over A G Tm and A: G {I,..., n} we get 

// ^l^P-mP dxdt = i // V" 1^(0 • Vl^P-mpd^dt, 

'J *7 Q'jp ^ J J Qt 

where 

Slkm = tt)), (t- 1A„T„M„<) , V< G R#^"*x("+!), 

the matrices T„, A„ and M„ being defined in §3.2.3. The lower triangular matrix 
A„ G Aln(®.) corresponding to the BDF2 method defined in Remark 3.7 satisfies 

(69) ||(A„)-i||, = ^(^I-^) <1 Vn>l, 

so that the time discretization fulfills Assumption (At). 

Finally, observe that thanks to the a priori estimates of §4.3 and in particular 
to Lemma 4.3, we find that there exists C depending only on uq, fl, T, q, d, and 
p* in (68) such that 

(70) [[ ?>;(,<p>dxdt<c||v(^p>||^,.(Q^), VpGcr(QT). 

J J Qx 

With the above “weak time derivative estimate” (70) and the “space derivative 
estimate” of Proposition 4.1 at hand, we can apply Proposition 3.8 and Theorem 3.9. 
We conclude that 


(71) 




n 

m 


U 


a.e. in Qt- 


m,n—^oo 
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Since n i® Uniformly bounded both in L°“((0, T); L^(n)) and in L'?+^((5 t) 

it follows from Proposition 4.1 and (66) that 

(72) {T^1nKn)m,n i® equi-integrable in L’'((0, T); L^(fl)) for all r G [1, oo). 

Applying Vitali’s convergence theorem we deduce the first claim of the following 
statement. 

Proposition 4.4. Let {Tm,£m,i^K)KeTm)m>i ® sequence of admissible dis¬ 
cretizations of n such that (68) holds. Let (n^)^ „ he the corresponding sequence 
of discrete solutions to the scheme (56), (58), and (59), then, up to an unlaheled 
subsequence, there exists u G L°“((0, T); L^(fl)) such that 

—> u strongly in L"^{{0, T)] L‘^{n)) for all r G [l,oo). 

m,n—¥oo 

Moreover, 4>{u) belongs to Lf{{Q,T)-,H^{M)) and 

V^{u) weakly tnL^Qrf. 

m,n—¥co 

The last claim of Proposition 4.4 is classical. Indeed, it follows from the first claim 
of the proposition that ) = (/)(7r” mJ)^) converge to </>(«); in addition, due 

to the uniform bound of Lemma 4.3, the discrete gradients converge 

weakly (up to a unlabelled subsequence) in L^{Qt) to some limit w. Then the 
weak limit w of f j can be identified with ’Vififu), by passing to the 

limit in the duality identities that express the action of V(/)(u) and <()(«(),) on 
test functions (see [27, 28] for details of this identification argument). Note that 
one can also prove that 

^ Vd\\vcf{u)\\mQ^y, 

prohibiting the strong convergence of </)(«]),) towards V^(m) if d > 2. 

4.5. Identification of the limit. The last step for proving the convergence of the 
scheme (56), (58), and (59) consists in passing to the limit in the appropriate weak 
formulation of the scheme, proving that the function u exhibited in Proposition 4.4 
is the (unique) weak solution to problem (50). 

Proposition 4.5. Let (frm,£m,i^K)Ke%n)m>i ® sequence of admissible dis¬ 
cretizations of LI such that (68) holds. Let ^ be the corresponding sequence 

of discrete solutions to the scheme (56), (58), and (59), then 

—>• u strongly in L''{{Q,T)-,Lf {LI)) for all r G [l,oo). 

m,n—^oc) 

where u is the unique solution to the problem (50). 

Proof. Let (p G C^{Lt x [0,T)) and let m,n > 1, then define = {pK)'^K^r^ 
by 

— -dx, VAT G 77,1, Vfc G {1,... ,n}. 

n J 

As in the proof of Lemma 3.10, we introduce the vector 

= {PK)°K^r^ = PIp. 


Pk = 


mK Jk 


P X, 
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the matrix A„ being defined by (32) and (34). It is straightforward to verify, thanks 
to the expression of and to the regularity of ip, that 

(73) Plv)\\L--{QTy —^ uniformly w.r.t. m. 

Multiplying the scheme (58) by and the scheme (59) by then 

summing over K S Tm and k € n}, reorganizing the sums, we find that 

(74) A^M+B^{p)+Cl{p) = 0, 


where 


Alip) = 




Qt 


^ TKLip'k-pl)). 

k=i KeT XLeAfK / 

C(¥^) =3 // V" V'(«;(.) • V" (^1 - pi,p) dxdt. 


Qt 

The discrete test-function has been built in order to ensure that A^{<p) can be 
rewritten as 


1 /" tt"’ P" 
m m m m 


pdxdt = - ff tt” m" 5” P” (/?dxd< -b [ 7r™M;(;°(x)v5(x,0)dx. 

j j (Dt J q 


' Qt Qt 

Therefore, using classical results (see e.g. [27]), we can easily check that 


(75) M(]j((p) —)■ -[[ udtpdxdt- f uop{-,0)d:>c, Vp € x [0,T)). 

m,n^oo 

It is now well known (see e.g. [30, 27] or [8]) that 

(76) B^^ip) [[ V^{u)-Vpdxdt. 

m.n^oo J Jq^ 

It follows from Lemma 4.3 and from (73) that 


(77) 


Clip) 


m,n—>00 


0 . 


Putting (75)-(77) in (74), we find that u is a weak solution of (50) Finally, as a direct 
by-product of the uniqueness of the limit value [41], one recovers the convergence 
of the whole sequence towards u in L‘^{Qt)- □ 
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